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With a help of the selfconsistent Hartree-Fock-Bogoliubov (HFB) approach with 
the DIS effective Gogny interaction and the Generator Coordinate Method (GCM) we 
incorporate the transverse collective vibrations to the one-dimensional model of the 
fission- barrier penetrability based on the traditional WKB method. The average fission 
barrier corresponding to the least-energy path in the two-dimensional potential energy 
landscape as function of quadrupole and octupole degrees of freedom is modified by the 
influence of the transverse collective vibrations along the nuclear path to fission. The set 
O , of transverse vibrational states built in the fission valley corresponding to a successively 

^ , increasing nuclear elongation produces the new energy barrier which is compared with 

C ' the least-energy barrier. These collective states are given as the eigensolutions of the 

GCM purely vibrational Hamiltonian. In addition, the influence of the collective inertia 
on the fission properties is displayed, and it turns out to be the decisive condition for 
, the possible transitions between different fission valleys. 

1. Introduction 

The phenomenon of nuclear fission is a complex process which involves the evolution 
of a nucleus over the multidimensional potential energy barrier emerged as a result 
of the competition between a short range attractive nucleon-nucleon interaction 
and the long range repulsive Coulomb forces. The problem of determining the most 
probable one-dimensional effective nuclear path to fission as an essential ingredient 
of the traditional WKB calculations of fission half-lives has been the subject of 
scientific efforts in low-energy nuclear physics since the time of the discovery of the 
nuclear fission phenomenon. 

It is clear that the commonly used WKB approach to estimate the fission barrier 
penetrability gives a discrepancy of fission half- lives up to a few orders of magnitude 
compared to experimental data. This fact points out on some essential deficiency of 
this approach. Above all, one should realize that the WKB approximation, based on 
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the assumption that the nuclear wave function evolves, starting from the ground- 
state or isomeric minimum, along a one-dimensional (least-action or least-energy) 
path toward the scission configuration in the multidimensional potential energy 
landscape, neglects the influence of transverse collective vibrations of the nuclear 
surface. These vibrations which are, in general, of different multipolarity, ampli- 
tude and frequency, depending on the actual nuclear configuration, always go along 
with the elongation of the nucleus on its way to fission due to the specific valley- 
like shape of the potential along the fission path. It is then obvious that the total 
collective energy is partitioned between all collective modes occurring along the nu- 
clear path to fission. Once the transverse collective modes are present, the effective 
fission barrier has to be understood as being determined by the set of the lowest 
vibrational states, built in the one-dimensional collective-potential well taken, for 
a given elongation, as the transverse cross section through the fission valley. 

A more fundamental approach would consist in assuming that the probability 
current related to the time-dependent nuclear wave function propagating essentially 
along the one-dimensional fission path has also non-zero components in all possible 
transverse directions. 

The above mentioned effect, by its influence on the fission-barrier shape, is then 
expected to modify the estimates of fission half-lives within the traditional WKB 
method in a substantial way. 

In the studies presented below we restrict ourselves to discuss the effect of trans- 
verse vibrations only along the least-energy asymmetric paths (with a substantial 
mass asymmetry of the fission fragments) which corresponds to the minimal po- 
tential energy as function of the elongation. This particular approach ensures that 
all vibrational states of the transverse mode possess higher energy than the corre- 
sponding collective-potential minimum. 

For a reasonable treatment of the collective properties of fissioning actinide 
nuclei the minimal set of collective variables has to comprise, at least, the nu- 
clear elongation, left-right asymmetry and possibly non-axiality. In the present 
work we study the least-energy paths, least-energy with transverse vibrations path 
and the least-action fission path for two actinide nuclei ^^^Th and ^ssy -^^j^jch in 
their ground states are practically stable against spontaneous fission. Already their 
shape-isomeric states can decay through spontaneous fission, a process which is 
competitive to the 7— back decay as shown in Rcf. ^. Since for those nuclei non- 
axial shapes are mostly important for configurations in the vicinity of the inner 
saddle point, (see e.g. Ref. the choice of only the elongation and the left-right 
asymmetry deformation degrees of freedom to uniquely describe nuclear states be- 
yond the isomeric minimum seems to be entirely justified. 

In order to take into account the collective modes the configurations mixing of 
pure Hartree-Fock-Bogoliubov (HFB) states in the total collective wave function 
is assumed. Applying the Generator Coordinate Method (GCM) together with the 
Gaussian Overlap Approximation (GOA) one determines the vibrational collective 
states by diagonalizing the GCM collective Hamiltonian. The knowledge of eigenso- 
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lutions of that Hamiltonian is essential for the problem of the transverse vibrations 
we are going to discuss in this paper. 

2. The Hartree-Fock-Bogoliubov method 

The potential energy landscapes and nuclear intrinsic states of the fissioning nu- 
cleus are determined in constrained HFB calculations with the DIS Gogny force ^ 
where the axially deformed nuclear states are described by quadrupole and octupole 
moments {92, 93} which create respectively the elongation and mass asymmetry of 
a nucleus. These intrinsic states J_*E'q) are defined in terms of quasi-particle vacuum 
states of the form (see e.g. Ref. ^ 

\^,) = l['iUo), 

i 

vt l*9> = 0, (1) 

where the quasi-particle creation operators rj^ are defined via the general Bogoli- 
ubov transformation. The explicit form of the operators r/^ can be deduced from 
the variational principle 

S{^{q2,q3)\H-XNN^XzZ-J2>^^Q^\^^iq2,q3))^0, (2) 

i 

where H is the many-body Hamiltonian built with the finite range effective Gogny 
force. The Lagrange multipliers Xn, Xz associated with the constraints on the 
neutron and proton number and the average nuclear deformations qi are calculated 
directly from the relations: 

{^{q2,q3)\N\-^{q2,q3)) = N, 
{^{q2,q3)\Z\'i'iq2,q3))=Z, 

(*(g2,g3)l4o|*(g2,g3)> = 

N+Z 

(Qio) = (^') = 0> (3) 

2 = 1 

where the last condition ensures that the center of mass of the system is fixed at 
the origin of the coordinate system. The multipole operators Q20 and Q^q in eq. ([3]) 
are defined as 



Q20 = a/ ^ > , rfY2o - > ^ {2zt -xf- yi) 
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(4) 



Equation ^ leads to the usual HFB equations which are solved by iterative 
diagonalization of the HFB Hamiltonian matrix expressed in the axial harmonic- 
oscillator double-centered basis until the required convergence is achieved. Con- 
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Fig. 2. Same as in Fig.[T]for the 238^ jm^leus. 
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sidering the symmetry properties of the HFB Hamiltonian one concludes that the 
entire potential-energy landscape is reflection symmetric with respect to (73 = 0. 
This directly implies that in both cases, ^^^Th and ^'^*U, presented in Figs. [1] and 
[2] one can immediately distinguish the symmetric and asymmetric potential energy 
valleys separated by the energy barrier appearing for well-elongated shapes. Both 
those valleys can be potentially treated as the nuclear fission paths leading either 
to the symmetric or (much more probable) asymmetric mass fragmentation (see 
Ref.lSfor 238u nucleus). 

In a purely static approach the criterion which determines the fission path is the 
potential energy value which has to be minimized over all deformation degrees of 
freedom entering eq. ^ except the leading parameter g2- Nevertheless, this static 
approach with the use of {(72, 93} as the deformation parameters can lead to artificial 
transitions between different possible paths if one is only interested in the absolute 
minimum of the potential energy with respect to q^. Clearly, (73 can vary only within 
the limited interval describing the pre-scission configurations shown in Figs. [U [2] 
with iso-energy lines. The absolute energy minima can occasionally belong to the 
symmetric valley, unlike the majority of elongated states as easily seen for several 
configurations of ^^^Th (i.e. from (72 ~ 190 6 to 92 ~ 230 b). Such sudden transitions 
between different fission valleys, which turn out to be highly probable in several Fm 
isotopes, strongly depend also on the collective inertia of the system and properties 
of the energy barrier between them. More reliable studies of such transitions can 
be performed with the help of the well known least-action formalism to be briefiy 
presented in section ID 



3. Beyond the mean-field approach 

We assume that the nuclear states are superpositions of the constrained HFB wave 
functions in the form given by the GCM theory (for review see 

I«'(92,g3)*> J j dq2dq3fi{q2,q3)\-^{q2,q3)), (5) 

where {92, 93} are the set of the collective variables identified as the quadrupole 
and octupole deformations respectively. The variational principle with a many- 
body Hamiltonian applied to the wave function defined by eq. ([5|) leads to the 
Hill- Wheeler equation (see Ref. ^ 

dq^dq'^ [H{q2, 93, 92: Id) ~ Ei /(g2, 93, 92) 93)] /j(92, 93) = (6) 
which yields the weight functions fi{q2, 93). Quantities entering eq. ^ are given as 

^^(92,93,92,93) = (*(92,93)|^|*(92,93)>, 

/(92, 93, 92, 93) = (*(92,93)|*(92,93)) (7) 

and denote respectively the GCM Hamiltonian and the norm kernel while Ei is the 
energy associated with the state |^'(92, 93)^). Instead of solving the integral equation 



9, 2008 0:48 WSPC/INSTRUCTION FILE kaz06a 



6 Authors' Names 

([6]) one can transform it to a much simpler form admitting the GOA approximation, 
(e.g. Ref. With this assumption eq. ([6]) reads 

'Hcoii9i{q2,qz) = Ecoiigt{q2,q3), (8) 

where now 54(92, 93) being the Gauss transform of fi{q2,q3) (see e.g. Ref. ^ is the 
collective wave function while Ti-coii is the collective Hamiltonian 

hj^ ^ d 1 d 

HcoU = V -^^—7 + VHFB('72,g3) - y] AVfe;(g2,g3) (9) 

fc^2 Bki[q2,q3) dqi 

which can be proved, in a straightforward calculations, to be formally identical 
to the collective Bohr Hamiltonian. The quantities T4jfb('J'2, 93), AVfci(g2,Q3), Bki 
are respectively the constrained HFB deformation energy, t he so-c alled zero-point 
(ZPE) vibrational and rotational energy corrections (see Refs.^l^^H|) and the inertia 
tensor associated with the quadrupole and octupole modes. The latter is calculated 

19 

in the Inglis cranking model . The ZPE vibrational correction origins directly 
from the fact that overlaps /, eq. ([7]) are not products of 5{q2 — q'2)S{q3 — q'^)- The 
rotational ZPE correction, in turn, is connected to the restoration on average of 
the rotational symmetry in the collective wave function \^{q2,q3)i). However, for 
the simplicity of the calculations the latter can be determined in an approximate 
manner as 



AKot = ^ ($(92,93) 



$('72,93)0, (10) 



2T 

n—1 ^ 

where j„ and Z„ denotes respectively the n*''— component of the total angular 
momentum and the nuclear moment of inertia with respect to the n*'' axis. 

The shortly presented GCM-f GOA procedure goes, obviously, beyond the stan- 
dard mean-field HFB approach providing in return, by the assumption of the con- 
figuration mixing, the average treatment of the long-amplitude nucleon-nucleon 
correlations which are not taken into account within the mean-field theory. 



4. Least-action approach 

The static approach presented in section [2] allowing for determining in an approx- 
imate way the nuclear path to fission is not entirely realistic since it neglects the 
collective inertia properties of the system. It was proved that the local behavior of 
the mass tensor can significantly change the fission path as compared to the one 
dete rmi ned by the energy-minimization condition. It is a well known fact (see e.g. 
Ref. that, in case the collective inertia is taken into account, the most favorable 
transition between two neighbouring configurations corresponding to slightly dif- 
ferent elongations passes through states for which either the energy or the inertia 
tensor decreases more rapidly. 
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Let us define the action in the two-d imensional potential energy landscape along 
the trajectory Li as e.g. in Refs. I ^'^ l l'^ l 

SiL,) - 2 J ^2[V{s) - E,oii]B{s)ds, (11) 

where s = s{q2, qs) describes parametrically a trajectory Li in the {52, 93} deforma- 
tion space. The points Si and S2 correspond to the turning points of this trajectory 
for which V{s) = EcoU- Within the GCM+GOA approach applied here the HFB po- 
tential energy is deepened by the zero-point vibrational and rotational corrections 
AVvib, AVrot respectively 

V = VHFB-AVrot-AV„i,- (12) 

It should be noticed that as far as the vibrational ZPE correction is concerned, it 
changes only insignificantly the shape of the fission barrier by shifting it down in 
energy by approximately the same energy for all minimum-energy configurations 
along the fission path, whereas the rotational correction lowers the barrier height 
(approximately) proportionally to the elongation. In the present calculations, it is 
assumed that the penetration of the fission barrier starts from the lowest vibra- 
tional state EcoU built in the isomeric potential-energy well and associated with the 
mode along the fission path. The projection of the inertia tensor Bij{q2, qs) on the 
trajectory Li is given in the standard way as 

B(^^T.'^.,i,.„>tt^ (13) 

ij=2,3 

The determination of the fission path corresponding to a more realistic treatment 
when both the energy and the inertia parameters are treated on the same footing 
would consist in finding the trajectory L^in which minimizes the action integral of 
eq. (fTT| . 

We performed the minimization of the action integral (jlip for two actinide nu- 
clei, ^^^Th and ^"^^U. As mentioned before, the resulting fission path is an interplay 
of the behaviors of the potential energy as well as the collective inertia. It is there- 
fore difficult, at first sight, to judge where the most probable nuclear path to fission 
passes through in the energy landscape studying only its local properties. How- 
ever, for both ^^^Th and ^38^ ^j^g symmetric and asymmetric energy valleys are 
relatively well pronounced (see Figs. [TJ [2|) the behavior of their least-action fission 
paths (LAPs) are remarkably different (see Fig.[3|). For the '^'^^Th isotope presented 
on the left hand side of Fig. [31 one observes that the LAP favors the configurations 
of substantially higher potential energies characterized by an about twice as much 
smaller inertia parameter (see r.h.s. panel of Fig. 0]) rather than the least-energy 
states. Certainly, in this nucleus the inertia effect plays a predominant role. In con- 
trast, in the right part of Fig. [3] we show a similar study for ^'^*U. The resulting 
LAP is now almost identical to the least-energy path. Obviously, for the latter case 
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Fig. 3. Potential energy barrier along the Lmin trajectory (thin) and the barrier corresponding to 
the least-energy trajectory (thick) for the actinide nuclei ^■^''Th and ^^^U. 




Fig. 4. (L.h.s. panel): the least-action trajectory Lmin (see text above) marked with thin line and 
the least-energy trajectory (thick line) obtained by minimizing the collective-potential energy V 
of eq. I I12I I over the octupole deformation 93, (r.h.s. panel): leading component of the mass tensor 
622 along Lmin trajectory (thin) and along the least-energy trajectory (thick) for ^■^^Th. 



the fission valley is deeper and stiffer as compared to the analogous asymmetric 
valley in ^^^Th. Hence, the static property of the potential-energy landscape, un- 
doubtedly, is of huger importance for this system, suppressing the influence of the 
collective inertia. 

The dynamic approach to the fission path can also be decisive for the problem of 
transitions between various energy valleys separated from each other by a relatively 
low energy barrier and corresponding to different fission modes. In fact, for the two 
actinide nuclei one observes only two possible valleys corresponding to the purely 
symmetric and asymmetric fission. Nevertheless, such a transition suggested for 
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•^^^Th by the static approach is finally excluded in the above sketched more realistic 
dynamical model. 

5. Collective vibrations across the fission path 

As mentioned in the introduction the nuclear fission as a collective mode can be 
strongly affected by a vibrational collective movement in the direction perpendicular 
to the fission path. 

The determination of the transverse mode corresponding to the symmetric fis- 
sion path consists in diagonalizing the one-dimensional collective Hamiltonian re- 
duced from the two-dimensional one, as expressed by eq. ([9]), by putting i = j = 3. 
For the majority of known actinide nuclei, however, a strong mass asymmetry of 
fission fragments is observed. The same treatment for an asymmetric path requires 
the knowledge of the local curvature of the path. On the other hand, since the 
constrained HFB calculations can be only performed on the discrete grid, the local 
curvature of the static fission path and the value of the inertia along the transverse 
path depend, to some extend, on the particular way of smoothing out the discrete 
energy values and the chosen mesh-point density. In order to wash out the discrete 
structure of the energy landscape as well as the three independent component s of 
the mass tensor we use the Gauss-Hermite folding method described in Ref. 
The necessity of applying an interpolation, naturally leads to an uncertainty of our 
estimates of the transverse-vibrations properties. 

Having determined the set of points corresponding to the least-energy fission 
path by minimizing the potential energy for a given parameter q2 with respect to 
53 within the asymmetric fission valley one determines, in turn, the transversal 
direction x which for any point on the fission path (given by {^2^13}) be 
locally expressed either as function of q2 or as function of 93 by the simple variable 
transformation 



where ip is the angle between the 53— axis and the direction. Knowing the deriva- 
tive of the fission path given as — f{q2) for given q2, we easily find the (^—values. 

The approach introduced here allows, by the particular choice of the two phys- 
ically relevant orthogonal collective degrees of freedom associated with the mode 
along the fission path and the transverse vibrations, for a more realistic prediction 
of the shape of the one-dimensional effective fission barriers which can serve as an 
input to the calculations of fission half-lives within the traditional WKB method. 
As a consequence of the separation of the mentioned two modes, the reduction of 
the full two-dimensional vibrational Hamiltonian ^ to the one-dimensional one of 
the transverse variable x and the derivative ^ is possible with the help of the trans- 
formation (|14p . In order to describe the vibrational excitations associated with the 




(93 - 93)^. 
(^2 - 92)51^, 



(14) 
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transverse collective movement, the reduced one-dimensional (GCM+GOA) collec- 
tive Hamiltonian reads 

where now = cos{f)-^, is the a;— component of the mass tensor which 
obeys the transformation rule given already in eq. (jl3p by replacing the curvilinear 
coordinate s along the barrier by the the transversal coordinate x 

= sm\ip) B22 - sin(2(^) B23 + cos'^{ip) B33. (16) 

The quantities VHFB(a;), ^V{x) are, as in eq. the HFB potential and the ZPE 
corrections taken along the x— direction. 
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Fig. 5. Static fission barrier obtained by the minimization of the potential energy given by eq. I)12|l 
with respect to 93 parameter (dashed Une), static barrier as a result of the potential-energy 
minimization along the perpendicular to the static fission path direction x for given q2— value 
(thick full line), the barrier with transverse vibrations (thin full line) 
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Figure [5] illustrates three fission barriers (all artificially shifted in energy to co- 
incide at the isomeric minima) obtained in different manners. The barrier marked 
with thin dashed line is obtained by minimizing the collective-potential energy 
^hfb((3'2, 93) + AV{q2,q3) (with AV denoting the sum of the vibrational and ro- 
tational zero-point corrections) with respect to the octupole deformation while 
the thick full line corresponds to the minimum of this energy as a function of 52 
measured along the direction perpendicular to the fission path. The third barrier 
represents the energy of the lowest vibrational states depending on (72 and resulting 
from the diagonalization of the one-dimensional collective Hamiltonian (|15p across 
the fission path. As one could expect, the last mentioned barrier including the ef- 
fect of transverse vibrations is effectively the lowest. This becomes clear when one 
realizes that the density of low-lying collective states built in the isomeric well is 
smaller as compared to the density in the saddle configuration. This directly implies 
that the distance between the collective potential energy minimum and the first vi- 
brational state of the transverse mode in the isomeric well is larger as compared to 
the corresponding distance for the outer-saddle configuration. 

Notice that the above mentioned effect can be easily explained in terms of the 
harmonic approximation where the energy interval between subsequent vibrational 
states is proportional to the vibrational frequency depending on the stiffness of the 
potential along the transverse direction in the vicinity of the isomeric minimum and 
the inverse of the effective mass parameter in this region. Therefore, the effective 
barrier height, measured as the difference between the lowest vibrational states of 
the isomeric minimum and the saddle configuration is, in fact, lower than for the 
static barriers shown in Fig. [5] by the order of 1 MeV for both investigated nuclei. 

6. Summary and conclusions 

We have presented three approaches to determine the effective fission path, 
demonstrating in this way a possible improvement of the commonly applied one- 
dimensional static approach to fission barriers. 

We have in particular confirmed in section [4] that, due to the balance between 
the static properties of the potential landscape and the collective inertia, significant 
changes can take place in the way nuclei to fission, particularly for nuclei charac- 
terized by relatively weakly pronounced energy valleys which determine the static 
fission paths on the average. The example of the nucleus ^^^Th has shown that in 
that case the least-action path leads through the states characterized by an about 
3-5 MeV higher potential energy than the least-energy path. Keeping in mind the 
one-dimensional WKB method, one easily concludes that the collective inertia is 
of the same importance for the barrier penetrability problem as the shape of the 
fission barrier. 

We have observed for the case of ^^®Th that a purely static treatment of the 
fission phenomenon would allow for transitions between different fission valleys, sep- 
arated by a substantially high energy barrier. In this particular case the rapid shift 
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between different fission modes would lead to two sudden modifications of the mass 
distribution of the fission fragments due to a slight variation of the elongation. Us- 
ing the method based on the action integral presented in section[4]one demonstrates 
that in ^^^Th such an unlikely transition does not exist. Nevertheless, according to 
our recent preliminary investigations, similar single transitions between asymmet- 
ric and symmetric modes might exist in the region of Fm isotopes. Investigations 
devoted to this topic are under progress. 

As already mentioned in the introduction, the one-dimensional description of 
the fission-barrier penetrability seems to be not sufficient since the theoretical pre- 
dictions of, for example, fission half-lives in different mass regions remain in dis- 
agreement with the experimental estimates by several orders of magnitude. Once 
the fission path is found we attempt in section 5 to extend this approach by con- 
sidering the influence of the collective transverse vibrations along this path using 
fully microscopic GCM-t-GOA approximation well suited for studies of heavy and 
super-heavy nuclei. Such an approach provides the consistent form of the collec- 
tive vibrational Hamiltonian without any adjustable parameters. As expected, the 
collective transverse movements cause a lowering of the effective one-dimensional 
fission barrier height by approximately 1 MeV for the discussed ^^^Th and ^ssy 
nuclei. This happens due to different static (potential energy stiffness) and dynam- 
ical (inertia) collective properties of the isomeric minimum and the saddle point of 
the fission valley. A collective effect of that order of magnitude can by no means be 
neglected in the estimates of fission half-lives or fission cross sections. 
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